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Abstract 



The analytic, nonlinear integral equation approach is used to calculate the finite-size 
corrections to the transfer matrix eigen-spectra of the critical dilute O(n) model on 
the square periodic lattice. The resulting bulk conformal weights extend previous 
exact results obtained in the honeycomb limit and include the negative spectral 
parameter regimes. The results give the operator content of the 19-vertex Izergin- 
Korepin model along with the conformal weights of the dilute Al face models in all 
four regimes. 



1 Introduction 

Among other physical phenomena, the integrable dilute O(n) model on the 
square lattice [1] is relevant to self-avoiding polymer chains in the bulk [2,3]. 
The partition function of the dilute 0(n) model is denned by [1,4] 



where the sum is over all configurations Q of non-intersecting closed loops 
which cover some (or none) of the lattice bonds. The possible loop configu- 
rations at a vertex are shown in Fig. 1, with a vertex of type i carrying a 
Boltzmann weight p». In configuration Q, rrii is the number of occurrences of 
the vertex of type i, while P is the total number of closed loops of fugacity n. 

The loop weights in (1.1) are [1,4] 
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Fig. 1. The 9 vertices of the dilute O(n) model. 
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Here n = — 2cos4A. These weights were originally constructed via a map- 
ping involving the Potts model [1] and later seen to satisfy the Yang-Baxter 
equation for loop models [4,5]. On the other hand, when mapped to a 3-state 
vertex model [1], the dilute O(n) model is seen to be related to the integrable 
19- vertex model of Izergin and Korepin [6]. The Nienhuis 0(n) model on the 
honeycomb lattice [7-9] follows from either of the special values u = A and 
u = 2A of the spectral parameter [1,10]. In the appropriate region the model 
thus contains the essential physics of the self-avoiding polymer problem at 
n = [2,7,8,11,12]. 

The dilute O(n) model has recently been used to construct a family of di- 
lute A-D-E lattice models [13-15]. These models are restricted solid-on-solid 
models with a finite number of heights built on the A-D-E Dynkin diagram. 
At criticality, the face weights are [13-15] 
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where the pi are as given above. The generalized Kronecker delta is unity if all 
its arguments take the same value and is zero otherwise. The Perron-Frobenius 
vectors S a in the face weights are the eigenvector of the largest eigenvalue of 
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the adjacency matrix A of the A-D-E graphs, 



Y,Aa,bS b = 2cos T —S a , (1.4) 

b L + 1 

where for the dilute Al models, L is the number of graph states, with a, b,c,d = 
1,2,- .-,L. 

The dilute 0(n) model exhibits various branches of critical behaviour [16-19]. 
These are reflected in the properties of the dilute A-D-E models, for which 
there are four physical branches [13] 
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Recent studies have highlighted the prominence of the dilute Al face models, 
which admit an off-critical extension [13,15]. In regime 2 the A% model lies in 
the same universality class as the Ising model in a magnetic field and gives 
the magnetic exponent 5 = 15 [13,15,20]. This A 3 model also shows the E 8 
scattering theory for massive excitations over the groundstate [21-23]. Both 
su(2) and su(3) fusion hierarchies of the dilute Al face models have been 
constructed in [24,25]. 

In this paper we both generalise and extend earlier calculations of the criti- 
cal properties, such as the central charges and bulk scaling dimensions (the 
conformal spectra), of the dilute O(n) model and the related dilute A L and 
Izergin-Korepin models. After outlining the necessary preliminaries in Section 
2, our calculations are presented in Section 3 for branches 1 and 2 and in 
Section 4 for branches 3 and 4. The method employed involves the extension 
of the nonlinear integral equation approach [26,27,18] to obtain the complete 
conformal spectra, as has been done for the six-vertex model [28,29] and most 
recently [30] for the Andrews-Baxter- Forrester (ABF) model [31]. Having read 
Section 2, those readers not specifically interested in the technical details may 
prefer to skip to Section 5 where a discussion of our results for the various 
models concludes the paper. 
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2 Bethe equations and known results 



As we are interested in bulk critical behaviour, we consider periodic boundary- 
conditions across a finite lattice of width N, where for convenience we take 
N even. The eigenvalues T(u) for the row-transfer matrix T(u) of the dilute 
0(n) model are given by [17-19] 



T{u) = e-'^ 



where 



s(2A-w)s(3A-w) Q(u + X) 
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(2.1) 



s(u) = sin N (u), Q(u) = Y\ cosh(iw — Uj) 

3=1 



(2.2) 



and the m zeros {uj} satisfy the Bethe equations 



ol<P 



cosh(-Uj + iA) 
cosh(-Uj — iA) 



yr smh(uj - u k + 2iA) smh(uj - u k - iA) 
^ sinh(-Uj — Uk — 2iA) sinh(-Uj — Uk + iA) 



for j — 1, . . . , m. It is convenient to label the sectors of T(u) by i = N — m, 
where £ = for the largest (groundstate) sector, £ = 1 for the next largest, 
etc. 

The Bethe equations ensure that the eigenvalues T(u) are analytic functions 
of u. Apart from the phase factors these equations are the Bethe equations 
of the Izergin-Korepin model [32-34]. In general is a continuous variable 
associated with a "seam" to ensure that loops which wrap round the cylinder 
carry the correct weight n. Thus 

= Ti - 4A (2.4) 



for the dilute O(n) model in the largest (£ = 0) sector of T(u) with <fi = in 
all other sectors. For the Izergin-Korepin model, = in all sectors. 

On the other hand, for the dilute A L face models there is a fixed number of 
Bethe roots (£ = 0) and 

= tts/(L + 1) (2.5) 
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with s = 1, . . . , L. In this case the transfer matrix T{u) has elements 



N 



(<7\T(u)\a>) = UW 









\<Tj 







(2.6) 



where the paths a = {<7i, 02, . . . , <Jn} and cr' = {a[, a' 2 , . . . , <j' n } are allowed 
configurations of heights along a row with periodic boundary conditions <Jn+i = 
U\ and cr' N+1 = a[. The face weights are those defined in (1.3). The expres- 
sion for the eigenvalues T(u) was given in [21] in terms of the more general 
elliptic functions. As we are also interested in the dilute O(n) model, we do 
not restrict the crossing parameter A to the values given in (1.5). This may 
lead to unphysical regimes in the dilute Al face model for which, however, the 
finite-size corrections to the transfer matrix eigenspectra are still of interest 
from the viewpoint of statistical mechanics and conformal field theory. 



2.1 Central charge 



Some exact results are known for the dilute O(n) model [17-19]. In particular, 
the central charge is found to be 



branches 1 & 2, (2.7) 
branches 3 & 4. (2.8) 

These results follow from the finite-size behaviour of the largest eigenvalue. 
The result (2.7) had already been obtained from the Bethe equations in the 
honeycomb limit [35-37]. However, the result (2.8) [18] had to await the de- 
velopment of the more sophisticated nonlinear integral equation approach 
[26,27,18] (see also [38,39]). 

The reason for this is that the distribution of Bethe roots for the largest 
eigenvalue differs significantly in each case. In the limit of infinite size N the 
Bethe roots are distributed on the lines [17-19] 



1 - 



30 2 



c = 



7T(7T-2A) 

3 3(f) 2 
2^A 



branches 1 and 2 3m (uj) = \n, (2.9) 
branches 3 and 4 Sm (uj) = ±|7rA. (2-10) 

Whereas there are no finite-size deviations from the line (2.9), the finite-size 
deviations from (2.10) are severe enough to render the more standard root 
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density approachF] invalid. Here we extend the analytic, nonlinear integral 
equation approach in the dilute O(n) model [18,19] to the calculation of the 
conformal weights in all four branches. Our treatment follows that given in 
the recent study of the ABF model [30]. 

The above results for the central charge have already been used to obtain the 
central charges of the dilute A—D—E models [13,15,19]. In particular, for the 
dilute Al face models, either (2.5) or the O(n) value (2.4) with the appropriate 
value of A in (1.5) gives 



1 - 



6 



h(h-l)' 
6 

h(h-l) 



branches 1 & 2, 



branches 3 & 4, 



(2.11) 



where 



L + 1, branches 2 & 4, 
L + 2, branches 1 & 3. 



(2.12) 



The first two branches give realisations of the unitary minimal series, while 
the other two branches involve a product of the unitary minimal series and an 
Ising model. 

The 0(n) model had earlier been identified [41,42] in the conformal classifica- 
tion scheme [43,44] with the aid of the Nienhuis Coulomb gas results [7,8]. In 
particular, c = 1 — Q(g — I) 2 /*?, where g G [1, 2], with g = h/(h — 1), in the high 
temperature phase (branch 1) and g G [0, 1], with g = {h — l)/h, in the low 
temperature phase (branch 2). Here g = 2(1 — 2X/tt). The Ising value c = ~ 
thus occurs both for the dilute A2 model (n — 1 in the high temperature 0(n) 
phase) and the dilute A 3 model (n = \/2 in the low temperature 0(n) phase). 
The central charges of the dilute A^ face models have recently been estimated 
numerically from the finite-size diagonalisation of the dilute A^ model trans- 
fer matrix for various L on all four branches [45] . The central charge has also 
been derived by solving the transfer matrix functional relations of the dilute 
Al model on branches 2 and 4 [46]. The calculation confirms the result (2.11) 
obtained via the dilute O(n) model [13,15,19]. 



2.2 Scaling dimensions 



Various scaling dimensions have been calculated via the Bethe equations for 
the dilute O(n) model. Again in the honeycomb limit for branches 1 and 2, 

1 See, for example, [40] and references therein. 
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the 'magnetic' set of scaling dimensions is found to be [35-37] 

X *~ 4vr(vr - 2A) " * 9 ' ~ " (2 ' 13) 



Alternatively, this result is written as 



2Ae/ 2 fl , branch 1 



X ° ~ ^ 2ATe'?> ! branch 2 



where 



A W _ - (ft - l)*] 2 - 1 

4/^-1) (2 - i5j 



is the Kac formula. This result had earlier been obtained via Coulomb gas 
calculations [47,48]. 

On the other hand, both numerical evidence [37] and root-density calculations 
[40] revealed the set of 'thermal' dimensions to be 

A -j(,-4A) = 2 . 
J 7T — 2A g 



Both the results (2.13) and (2.16) generalized earlier results via the Coulomb 
gas [7, 8, 11]. 2 The thermal dimensions follow from 

ye _ / 2A i,2j+i , branch 1 
j \2A 2i+lil , branch 2 1 U) 



in the Kac formula [37]. 

On the other hand, the situation is not so clear for branches 3 and 4 of the di- 
lute square lattice model. Numerical evidence [18] indicates that the magnetic 
dimensions are given by 

t 2vr 8ttA V ' 



The only known thermal result is X{ = 1 [18]. There are no Coulomb gas 
results for these branches. 



2 The leading thermal dimension had been conjectured earlier for the O(n) model 
by Cardy and Hamber [49]. 
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The conformal weights of the dilute Al face models have been estimated nu- 
merically from the finite-size diagonalisation of the transfer matrix (for L = 3 
and L = 4 at u = 3 A/2 on all four branches) [45] and from numerical solutions 
to the Bethe equations for L = 3 [23]. For branches 1 and 2, the results fulfil 
the expectation that the scaling dimensions reflect the conformal weights of 
the unitary minimal series. For branches 3 and 4, they reflect a product of 
the Ising and unitary minimal series. The related modular invariant partition 
function has been discussed at length in [45]. 

As mentioned above, the finite-size corrections to the transfer matrix eigen- 
spectra have been obtained for the dilute face models in branches 2 and 
4 in [46] via the functional relation method [38,39,50]. The analytic calcula- 
tion confirms the conformal weights obtained via the calculation of the local 
height probabilities for L odd [20]. Here we consider the dilute models in 
all four branches with more general crossing parameter A and calculate the 
conformal spectra for each branch. 



3 Branches 1 and 2 



We consider branches 1 and 2 defined by 

0<m<3A ti/6 < A < 7r/3. (3.1) 

This regime covers the A values (1.5) for the dilute A L models. However, the 
derivation below is also valid for the dilute O(n) and Izergin-Korepin models 
in the larger interval < A < 7r/3. Let us introduce the new variable v — iu 
with a shift Vj = Uj — |i7r. The Bethe equations (2.3) are then of the form 

P(vj) = -1, (3.2) 

where 

P[V) ~ 6 $(v + i\)q(v + i\)q(v - 2iA) ' ^ 

m 

®(y) —smh N v , q(v) = J[ smh(v - Vj) . (3.4) 

After the shift, the Bethe roots Vj are distributed along the real axis, with 

q(v) = q(v) and p(v) = l/p(v) (3.5) 

where the overbar denotes complex conjugation. 
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3.1 Nonlinear integral equation 



Define two functions that are Analytic and Non-Zero (anz) in the strips 
around the real axis: 



a(v) = e lu g(v) p(v + i£), 

A(v) = l + a(v)/g(v). (3.6) 

The phase factor u has been introduced for taking different branches of the 
log function involved in the subsequent Fourier transforms. We take 



u = < 



' sgn(v)(£ — r)n dilute O(n) model, t ^ 
rvr dilute O(n) model, i = (3.7) 

k 7r(r — s) dilute A L face model 



where the function sgn(t> ) = —1 for Re(v) > and +1 otherwise. For the 0(n) 
model the integers r, s are restricted (as discussed further in section 5). For 
the moment we leave them as arbitrary integers. The function g is introduced 
for compensating the anticipated bulk behaviour of p(v + i£) and is given by 

g(t,) = V tMT,-iA + ifl J' P = 7F/ ' ( } 



where < ^ < \n. The i7T— periodic function a can be rewritten as 

= ^ } *(«■ - U + mv - iA + |+ ■ ifM, + 2.A + if) 
v ; yv J ®(v + i\ + i£)q(v + i\ + i£)q(v-2i\ + m + i£) v ; 



The above treatment results in the function a representing finite-size correc- 
tions. To see this we consider the Fourier transform pair 



1 If \" ■ 

a(k) = — J (\na(v)j e~ lkv dv 



— oo 

oo 

lna(u)) = / a{k)e ikv dk (3.10) 



for a and similarly for A. The Fourier transform of q(v) is defined to be 



oo+lr 



q(k) 



2tt J 



In q(v) 



-ikv 



dv < r < n, 



— oo+ir 
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n oo 

(lng(u)) = J q(k)e [kv dk < 3m (v) < ir. (3.11) 

— oo 

To represent a(k) by A(k) and A(k) we also need another relation, which can 
be given by applying Cauchy's theorem to the auxiliary function 

h(v) = i±PM, (3.12) 
p(v)q(v) 

which satisfies the non-trivial analyticity property 

oo+i£ oo— 

1 (ln/i(u)) e~ ifc,; civ = y (m e - ito rft;. (3.13) 



-oo+i£ — oo— 



From the equations following Fourier transforming (3.9) and inserting (3.12) 
into (3.13) we obtain 



Nke { z nk) cosh{\\k) 
g( )_ 2smh(±7rfc)cosh(§Afc) 

+ ^P^^^i^m ~ e-*m\ (3.14) 

2sinh(i7rA; - AA;)cosh(|AA;) V V; V V K J 

a{k) = F{k) A(k) - F^k)A(k), (3.15) 



where 



Fc(k ) - 2gfc sinh(A£;) cosh(|7r£; - |A/c) 

WJ ~ cosh(|AA;)sinh(|7rA;-AA;) 1 J 



and F(A;) = F (k). Transforming back and integrating twice we obtain the 
nonlinear integral equation 

\na(v) = F *\nA- F^lnA + C + C'v, (3.17) 
where the convolution is defined by 

oo 

(f*g)(v)= J f(w)g(v-w)dw. (3.18) 



The constant C is chosen to be C = for all terms to remain finite. The 
other constant C is scaling-dependent and is fixed after taking the scaling 
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limit defined by 



a±(x)= lim a(±v)/g(±v), 

N— »oo 

A±(x)= lim A(±v) = l + a±(x). (3.19) 

AT— >oo 

The nonlinear integral equation then becomes 



in which the kernel if is given by 

(3 ' 21) 



where 



=7 2 (-n-) = ^(±-^ + 2i{). (3.22) 

We can see that K T (x) = K(—x), a key property to be used in the derivation 
of the finite-size corrections. Taking x — > oo we obtain 

C± = i7r(w±-^)/(7r-2A), (3.23) 



where 



f c<j =p 2£A dilute O(n) model, ,„ „a 

^± — 4 ti / r , i (3.24) 

La; dilute A L face model 



The nonlinear integral equation (3.20) is equivalent to the Bethe equations for 
the largest eigenvalue, as given in [18,19]. The key difference is the change in 
the integration constants C± for the low-lying excited states. This is similar 
to the nonlinear integral equation approach in the six-vertex [28] and ABF 
[30] models. In each case the constants contain the necessary information to 
extract the conformal weights. 



3.2 Conformal spectra 



The eigenvalues of the transfer matrix are given by (2.1). For small positive 
values of u the first term in the eigenvalue expression dominates exponentially. 
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For small positive imaginary v we therefore have 

T{v) ~ e-^(v - 2i\)<S>(v - 3iA) ^^|^ (3.25) 



for the finite-size corrections. Taking Fourier transforms and integrating twice 
yields 



, r^, n ~\ t f / \ 2\/3p °f /sinh4p(t> — w — if) , 

111^=-^,;+— / -wTT ^A(w 

7r J \^smh6p(i> — w — i£J 



sinh 6p(f — w + i£) 
where the free energy is given by 



sinh4p(t> — w + if) , _ \ 

^ ; \nA(w) dw, (3.26) 



/• sinh(/cif ) sinh(3/cA + if k) cosh(5A;A — kit) cosh(A;A) 
f ^ v) = 2 J M A;cosh(3A;A)sinhfM ^ 



-co 



The integration constants have been fixed again by the limit v — > oo. Taking 
the thermodynamic limit N — > oo in (3.26) and using the definitions (3.19) 
gives 



\nT(v) = -Nf oc (v) + ^^e 2 ^m (e' 2 ^ J \nA + {x)e~ x 

\ — oo 

up to order 1/iV. To proceed further we consider the expression 



oo / 



dx (3.29) 



which can be written exactly as 



L(z) + L(l/z) = ^- (3.30) 
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in terms of the Rogers dilogarithmic function 



L(x) 



/ln(l +y) Iny 



\ y i + y 



dy. 



(3.31) 



On the other hand, making use of (3.20) in (3.29) and using a±(— oo) 
a±(— oo) = and a±(oo) = l/a±(oo) = e 1 ^ - ^, we arrive at the result 



+ 8v x 3.w/ [ < : 



7r(^± 



7T-2A 



(3.32) 



Equating (3.32) and (3.30) gives the integral in (3.28). Thus inserting this 
integral into the expression (3.28) we obtain 



, r^, n t* y /■ / \ vr sin(2ipt>) . ^, AS 
ln7» = -7V/o» ^^(c-24A) 



(3.33) 



to leading order in 1/N. This is our final result, from which the central charge 
and conformal weights can be read-off [51-53] as 



c=l - 



30 2 



A= <^ 



7r(7T-2A)' 

f (w - =F 2A) 2 - (tt - 4A) 2 
87r(7r - 2A) 
- 0) 2 - (tt - 4A) 2 
87r(7r - 2A) 



dilute O(n) model, 
dilute At face model. 



(3.34) 
(3.35) 



4 Branches 3 and 4 



On branches 3 and 4 the spectral and crossing parameters are specialized in 
the regime 



-7T + 3A<M<0 



1 , 1 

-71 < A < -7T 

6 ~ 3 



(4.1) 



The following computation of the finite-size corrections to the transfer matrix 
eigenspectra for each of the models is valid for the larger interval < A < \u. 
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We proceed in a similar manner as for branches 1 and 2 and introduce a new 
parameter v = iu and set Vj = Uj. The function p(v) is defined by 



AJ $(t; + iA + ±7ri)g(t; + iA)g(t;-2iA)' 



(4.2) 



with $ and g as given in (3.4). In [18,19] it has been checked that the Bethe 
ansatz roots are (almost) located on the lines 3m (v) = ±|A in the complex 
v-plane. As a consequence we still have the symmetries of equation (3.5). 



4-1 Nonlinear integral equation 

We proceed by defining functions that are anz in the strips around the real 
axis: 

A{v) = 1 + a(v)/g(v) a(v) = g(v)p(v — iA) [1 + p(v)} 

p(v)p(v — iA) 



B(v) = l + 0(v)/g(v) P(v)=9W 1 + l){r _ jA) 

C{v) = 1 + l(v)/g{v) 7(u) = g(v)p(v - iA) 

5(f) = p(v). 



(4.3) 



The function g{v) = th p(t> + iA — |vri), with p = 7r/ (27r — 6A), is introduced 
to compensate the anticipated bulk behaviour of the functions a,/3, 7. 

We define the Fourier transform of the functions a, (3, 7 as in (3.10). For q we 
have 



1 


oo+ir 


27 


/ 




-oo+ir 


(in 


q(v)j 


1 


oo+ir 


27 


/ 



-ikv 



-oo+ir 



\nq(v)j = J qi (k)e ikv dk 



-7T + \l\ < r < — |7T 



g(fc)e ito rfA; -vr + ±tt < 3m (u) < — |tt (4.4) 



-±7T < r < |?r 



■|7T < 3m (t>) < |7T. 



(4.5) 



To solve the functional relations we need the relations of the Fourier transforms 
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of a, (3, 7. First we can see that not all functions are independent and thus we 
have 



(3{k)- 1 {k)-5{k) + C{k) = 
a(k) - a(k) - 7 (Jfe) + j(k) - 5(k) = 
A(k) - B(k) - C(k) = 0. 

Applying the Fourier transform to the 5, 7 gives 



(4.6) 



7(fc) = iV/c sinh A/c/ sinh vr£;/2 - \Nke^ / cosh -(3A - vr) 



+ (e- 



-Afc+7rfe 



+ e 2Afc -e 3A )g(A;)-g 1 (£ ; ) 



(4.7) 



k k 
S(k) = Nk sinh Xk/ sinh 7rfc/2 - 4e^ sinh Afc/2 cosh -(3A - vr) q(k). (4.8) 

Other relations follow by applying Cauchy's theorem to the auxiliary functions 



h 1 (v)=p(v- ±iA)[l+p(u + ±iA)], 

l+p(t,-liA)[l+p( V + ^iA)] 

2 W — 7 — ; i-u ) 

p(v + 21A) 

/i 3 (u) = [1 +P(^ - ±i\)[l+p(v + ±iA)]/g(u - §iA), 
which all satisfy the non-trivial analyticity property 



(4.9) 
(4.10) 
(4.11) 



\nq(v)) e~ ikv dv 



00— ^lA 

/ (' 

-00— ^lA 



In q(v) ) e ikv dv . 



(4.12) 



It follows, respectively, that 



a(k) = -e Xk (3(k), 
A{k)-5{k) = e Xk [A{k) + 5{k)], 



(4.13) 
(4.14) 



and 



e-2 Xk (A(k) - e Xk q(k)) 
(_ 

B(k) - (3(k) ~ Qi(k) + 



e2 



A A; 



V 



\Nke~2 Xk 
cosh(|AA; — \nk) 



(4.15) 
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Now solving (4.6)-(4.15) and their complex conjugates in terms of the func- 
tions A and B, we find 



a(k) + 7 (Jfe) = F(k)A(k) + G(k)A(k) + H(k)B(k) + H{k)B{k) 
(3(k) - 7 (Jfe) = H(k)A(k) + H(k)A(k) + B(k) 

C(k) = A(k) - B(k) (4.16) 



Nke 2 ^ cosh 



q(k) 



Xk 



e~2 {n+x)k A(k) - e-2 (7T ~ x)k A(k) 



with 



F(k) 
G{k) 
H{k) 



2sinh ^7r/ccosh |(3A - n) 4 sinh A/c cosh §(3A - 7r) 



sinh |(tt - 3A) - 2 sinh |(tt - 5A) 

2 sinh AA; cosh |(3A - 7r) 
2 e -|(7r-5A)fe _ 2 e -|( 7r - 7A ) fe _ 2 e _ |( 7r_3A ) fc -|- e |( 7r - A ) fe 
4 sinh A A; cosh | (3 A - 7r) 



e 2 



AA- 



2 cosh | A A; 



(4.17) 



(4.18) 



Transforming back and integrating twice, we obtain a coupled set of nonlinear 
integral equations, 

\na(v) + ln 7 (-y) = F*\nA + G*\nA~ + H*lnB + II*lnB~ + C, 
lnP(v) -\n-f(v) = H *\nA + H *\nA + \nB, (4.19) 

where we have introduced 



F(v) = — I F(k)e lkv dk, 
2n ' 



— oo 
oo 



G(v) = ^- I G{k)e lkv dk. 



(4.20) 



-oo 
oo 



H{v) = [ H(k)e ikv dk. 

— oo 

Taking the scaling limit as in (3.19), the nonlinear integral equations become 



/ In a± + In c± \ 
In b± — In c± 
In a± + In c± 

V \nb± — lnc± / 



±4ie~ 





_ e TpiA 

V o J 



+ K * 



/In At \ 
ln£± 
\nA± 

\hxB ± J 



( 1 


-1 
V o 



,(4.21) 
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where the kernel K again satisfies the very useful symmetry K T (x) = K(x). 
The integration constant in (4.19) follows from the limit x — > oo. We have 

C± = iri(uj ± - <j>)/(2\). (4.22) 



Again the integration constants C± contain the essential information to obtain 
the conformal spectra. 



4-2 Conformal spectra 



For branches 3 and 4 the transfer matrix eigenvalues in (2.1) are dominated 
by 



T(v + iX- ivri) ~ j*Q(v + iA + ±7ri)$(<; + Wi) ^ ^ A(v). (4.23) 
Taking Fourier transforms and using the solution (4.16) for q(k), we obtain 

InT^ + iA-Iyri) = -NfM 



00 



+ £ / ( -J^W _ -J^W _ ) A,, (4.24) 

7r j ycosh2p(f — w + |iA) cosh2p(t> — w — p\) J 



-00 



where the bulk free energy is given by 



1 \ 9 f au srnn (^i' 1 ') sinh(3A;A + ivk — nk) cosh(5A;A — kir) cosh(A;A) . . 

Taking the thermodynamic limit iV — > 00 and using the definition as in (3.19) 
gives 



lnT(u) = -7V/oo(u) + j^e 2(m $fm e p[x J In A^e - * 

V -00 y 

+-^e- 2 ^9fm ^e-" iA J In A.^e^j dx. (4.26) 
To calculate the integral, we consider the expression 
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/ In a± + In c± \ ' / In A± \ T / In a± + In c± \ / In A± \ 



In — In c± 
ln a± + In c± 



lnB± 
lnXt 



\ln6±-lnc±/ Vln£ ± / 



In b± — In c± 
In a± + In c± 



lnB± 
lnXt 



\ln&± -lnc±/ \lnfi±/ 



rfx,(4.27) 



which can be evaluated exactly using the Rogers dilogarithmic function rela- 
tion (3.30). We thus arrive at 



L(a±(oo)) +L(Va±(oo)) + l(i/6±(oo)) 

+L(b±(oo)j + L^c±(oo)) + l(i/c±(oo)) = vr 2 - 8A;tt 2 , (4.28) 

where = 0, 1 [28]. We have also used the asymptotics of the functions 

a±(oo) = e i ( w± -*)(e i ( w± -^ + 1), 6±(oo) = e 2 ^±"^/(e^±"^ + 1), c±(oo) = 
e i(w±-«) and a± (_oo) = c±(-oo) = c±(-oo) = 0. 

On the other hand, substituting (4.21) into (4.27) we arrive at 



163m (e ±pix J \nA±(x)e- x ) + 



■k(uj± - 4>)' 
A 



(4.29) 



Combining the results (4.28) and (4.29) we are left with 



3m \e ±plX 



J hi A±(x)t 



Inserting this in the expression lnT(f) we obtain the final result 



ln7> + iA - ivri) = -NfM + ^^(c - 24A). 
The central charges and conformal weights are given by 



(4.31) 



3 3(vr-4A) i 



2 2ttA 

( (cj-0^2£A) 2 - (tt-4A) 2 



A=<^ 



16ttA 

(u - 0) 2 - (tt - 4A) 2 
16ttA 



+ Ai S i n g dilute O(n) model 



(4.32) 



(4.33) 



+ A 



Ising 



dilute Al face model 



with A Ising = 0, \. 
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5 Summary and discussion 



We have calculated the finite-size corrections to the transfer matrix eigen- 
spectra of the intimately related dilute 0(n), dilute Al and Izergin-Korepin 
models at criticality via the nonlinear integral equation approach. The re- 
sulting conformal weights defining the critical exponents are seen to follow 
from appropriate branches of the log functions appearing in the Bethe equa- 
tions. For the dilute O(n) model the integration constants appearing in the 
key nonlinear integral equations (3.20) and (4.21) differ in the distinct limits 
v — > ±oo, leading to two constants, C±, in the scaling limit. However, they 
satisfy |C+| = |C_|, allowing the calculation to go through in a similar manner 
as for the ABF model [30]. 



5.1 Dilute 0(n) model 



Consider first the dilute O(n) model in branches 1 and 2 in the u positive 
regime. Our final results are the central charge (3.34) and the conformal 
weights (3.35). These are in agreement with the previous results outlined in 
Sections 2.1 and 2.2 with the 0(n) <p value (2.4). In particular, the magnetic 
dimensions (2.13) follow with the parameters £ ^ and r = in (3.24). Simi- 
larly, the thermal dimensions (2.16) follow with £ = and r = 2j. As expected, 
in this regime the conformal dimensions are seen to be in agreement with the 
results obtained in the honeycomb limit. 

On the other hand, in branches 3 and 4 of the u negative regime our final 
results are (4.32) and (4.33). Here the conformal dimensions are new. The 
conjectured magnetic dimensions (2.18) [18] are associated with £ ^ and 
r — t. The similarly conjectural thermal dimension X{ = 2Ai S i ng = 1 follows 
from the choice £ = and r = 0. More generally we see that this dimension 
belongs to the thermal set 



_ J 2 7T-i(7T-4A) 

A j+1 - Rising, (0-34) 

which is given by setting £ = and r = 2j in (4.33). 



5.2 Dilute Al face model 



Recall that the Bethe equations of the dilute A L face model follow from the 
choice of crossing parameter given in (1.5) with seam <fi as given in (2.5). In 



19 



this case the branches 1 and 2 results (3.34) and (3.35) give the central charge 
(2.11) and conformal weights (2.15) of the unitary minimal series. 

In a similar manner, the branches 3 and 4 results are indicative of the product 
of the unitary minimal series with the Ising model. The central charge is given 
by (2.11) and the conformal weights again by given by (2.15), with however, the 
additional A Ising component. Here A Ising = 0, |, ^ is expected, in accordance 
with the results of [20]. However, ^ does not appear in our results. Our method 
needs further refinement to reveal this conformal dimension. 



For the more general crossing parameter 

A = - 1 + 



4 V L + Y 

with \k\ < [(L + 1)/3J, the integer part of the fraction (L + l)/3, our results 
(3.34) and (4.32), and (3.35) and (4.33) imply 



6k 2 

c — 1 — — — — branches 1 and 2 (5.35) 

(L + 1)(L + 1 - k) 

3 Qk 2 

c— — r— — branches 3 and 4 (5.36) 

2 (L + l)(L + l-k) v ; 

for the central charges and 



k 2 

— branches 1 and 2 (5.37) 

k 2 

h Ai s i ng branches 3 and 4 (5.38) 

for the conformal weights, where s — 1, 2, • • • , L and t — 1, 2, • • • , L — k. These 
results indicate the non-unitary minimal models for the dilute Al models in 
braches 1 and 2 and the non-unitary minimal models plus an Ising model in 
branches 3 and 4. For the case k — 1 our results confirm the conformal weights 
presented in [20]. For k > 1 our results show that the dilute A L models in 
branches 1 and 2 can be classified by the same universality classes as for the 
ABF models [54,30]. 



[(L + l)t- (L + l -k)sf- 
4(L + l)(L + l-fc) 

l(L + l)t-(L + l-k)s] 2 - 
4(L + l)(L + l-fc) 



5.3 Izergin-Korepin model 



Our results for the O(n) model reduce to those of the Izergin-Korepin model 
when the seam = 0. In this way the central charge and conformal weights 
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are given by 



£ 2 (vr-2A) r 2 vr 



for < u < 3A with < A < 7r/3. On the other hand, 



3 v A£ 2 m 2 oA 
c=-, X e>m = — + — + 2A Ising (5.40) 



for — 7r + 3A < u < with < A < 7r/3. The result (5.39) is in agreement with 
the previous results in the honeycomb limit, while the result (5.40), reflecting 
also the additional Ising content, is new. 



It is a pleasure to thank Ole Warnaar for helpful discussions over the years. 
This work has been supported by the Australian Research Council. YKZ also 
thanks the Natural Science Foundation of China for partial support. 



References 



[I] B. Nienhuis, Int. J. Mod. Phys. B 4 (1990) 929. 

[2] B. Nienhuis, in Fundamental problems in statistical mechanics VII, H. van 
Beijeren ed. (North-Holland, 1990) pp 255-271. 

[3] B. Nienhuis, Physica A 163 (1990) 152. 

[4] S.O. Warnaar and B. Nienhuis, J. Phys. A 26 (1993) 2301. 

[5] B. Nienhuis, S.O. Warnaar and H.W.J. Blote, J. Phys. A 26 (1993) 477. 

[6] A.G. Izergin and V.E. Korepin, Commun. Math. Phys. 79 (1981) 303. 

[7] B. Nienhuis, Phys. Rev. Lett. 49 (1982) 1062. 

[8] B. Nienhuis, J. Stat. Phys. 34 (1984) 731. 

[9] R.J. Baxter, J. Phys. A 19 (1986) 2821. 

[10] N.Yu. Reshetikhin, J. Phys. A 24 (1991) 2387. 

[II] B. Nienhuis, in Phase Transitions and Critical Phenomena, Vol 11, C. Domb 
and J.L. Leobowitz eds. (Academic, 1987) pp 1-53. 

[12] B. Duplantier, in Fundamental problems in statistical mechanics VII, H. van 
Beijeren ed. (North-Holland, 1990) pp 171-223. 

[13] S.O. Warnaar, B. Nienhuis and K.A. Seaton, Phys. Rev. Lett. 69 (1992) 710. 



21 



[14] Ph. Roche, Phys. Lett. B 285 (1992) 49. 

[15] S.O. Warnaar, B. Nienhuis and K.A. Seaton, Int. J. Mod. Phys. B 7 (1993) 
3727. 

[16] H.W.J. Blote and B. Nienhuis, J. Phys. A 22 (1989) 1415. 

[17] M.T. Batchelor, B. Nienhuis and S.O. Warnaar, Phys. Rev. Lett. 62 (1989) 
2425. 

[18] S.O. Warnaar, M.T. Batchelor and B. Nienhuis, J. Phys. A 25 (1992) 3077. 

[19] S.O. Warnaar, Solvable loop models and graphs, PhD thesis, Amsterdam (1993). 

[20] S.O. Warnaar, P.A. Pearce, K.A. Seaton and B. Nienhuis, J. Stat. Phys. 74 
(1994) 469. 

[21] V.V. Bazhanov, B. Nienhuis and S.O. Warnaar, Phys. Lett. B 322 (1994) 198. 

[22] S.O. Warnaar and P.A. Pearce, J. Phys. A 27 (1994) L891. 

[23] U. Grimm and B. Nienhuis, Scaling properties of the Ising model in a field, 
|hep-th/9511174| . 

[24] Y.K. Zhou, su(2) fusion hierarchies of dilute models, Int. J. Mod. Phys. B (in 
press) . 

[25] Y.K. Zhou, P.A. Pearce and U. Grimm, Physica A 222 (1995) 261. 

[26] A. Kliimper and M.T. Batchelor J. Phys. A 23 (1990) L189. 

[27] A. Kliimper, M.T. Batchelor and P.A. Pearce, J. Phys. A 24 (1991) 311. 

[28] A. Kliimper, T. Wehner and J Zittartz, J. Phys. A 26 (1993) 2815. 

[29] M. Fujimoto, J. Phys. A 28 (1995) 1519. 

[30] Y.K. Zhou, J. Phys. A 28 (1995) 4339. 

[31] G.E. Andrews, R.J. Baxter and P.J. Forrester, J. Stat. Phys. 35 (1984) 193. 

[32] N.Yu. Reshetikhin, Sov. Phys. JETP 57 (1983) 691. 

[33] V.I. Vichirko and N.Yu. Reshetikhin, Theor. Math. Phys. 56 (1983) 805. 

[34] V.O. Tarasov, Theor. Math. Phys. 76 (1988) 793. 

[35] M.T. Batchelor and H.W.J. Blote, Phys. Rev. Lett. 61 (1988) 138. 

[36] J. Suzuki, J. Phys. Soc. Japan 57 (1988) 2966. 

[37] M.T. Batchelor and H.W.J. Blote, Phys. Rev. B 39 (1989) 2391. 

[38] P.A. Pearce and A. Kliimper, Phys. Rev. Lett. 66 (1991) 974. 

[39] A. Kliimper and P.A. Pearce, J. Stat. Phys. 64 (1991) 13. 



22 



[40] J. Suzuki, T. Nagao and M. Wadati, Int. J. Mod. Phys. B 6 (1992) 1119. 

[41] Vl.S. Dotsenko and V.A. Fateev, Nucl. Phys. B 240 (1984) 312. 

[42] V. Singh and B.S. Shastry, Pramana - J. Phys. 25 (1985) 519. 

[43] A. A. Belavin, A.M. Polyakov and A.B. Zamolodchikov, J. Stat. Phys. 34 (1984) 
763; Nucl. Phys. B 241 (1984) 333. 

[44] D. Friedan, Z. Qiu and S.H. Shenker, Phys. Rev. Lett. 52 (1984) 1575. 

[45] D.L. O'Brien and PA. Pearce, J. Phys. A 28 (1995) 4891. 

[46] Y.K. Zhou and P. A. Pearce, Conformal weights of dilute A lattice models (1995) . 

[47] H. Saleur, J. Phys. A 19 (1986) L807. 

[48] B. Duplantier, J. Stat. Phys. 49 (1987) 411. 

[49] J.L. Cardy and H.W. Hamber, Phys. Rev. Lett. 45 (1980) 499. 

[50] A. Kliimper and P. A. Pearce, Physica A 183 (1992) 304. 

[51] H. W. J. Blote, J. L. Cardy and M. P. Nightingale, Phys. Rev. Lett. 56 (1986) 
742. 

[52] I. Affleck, Phys. Rev. Lett. 56 (1986) 746. 
[53] J.L. Cardy, Nucl. Phys. B 270 (1986) 186. 
[54] P. J. Forrester and R. J. Baxter, J. Stat. Phys. 38 (1985) 435. 



23 



